In the early universe chemical equilibrium between particles like gluons and 
quarks was sustained by annihilations. On dimensional grounds, the 2 <-> 2 
annihilation cross sections at energies much higher than the masses of the 
particles involved should decrease as g 4 (fx)s~ 1 ~ g 4 (T)T~ 2 , where g(fi) is the 
coupling renormalized at the scale \i. Thus at very high temperatures, and 
in a radiation dominated universe, the total averaged annihilation rate should 
scale as 1"*°^ ~ g 4 (T)T. One should then compare r*nn with the Hubble rate 
H = l/2t ~ T 2 , and it is rather obvious that at sufficiently high temperatures, and 
for sufficiently weak interactions, 2^2 annihilations could not have maintained 
<^ \ equilibrium. For instance, given only the QCD interactions which become, on the 
average [1], weaker as T — > oo, at some high temperature quarks and gluons may 
not have been in thermal contact at all. If so, it certainly would have interesting 



ramifications for our view of the very early universe as ideal gas in thermal 
equilibrium. In fact, it was already roughly estimated in [2] that in SU(5) grand 
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unified theory equilibration is possible only at temperatures T < 3 x 10 15 GeV, and 
if the initial state is far out of equilibrium, chemical equilibrium is established 
at temperatures much lower than this [3]. Whether a grand unified theory 
actually exists or not is of course unknown, but it is very likely that QCD 



interactions were in full thermal equilibrium prior to the QCD phase transition 
at T = Tqcd — 200 MeV. Given that, an interesting question then is, within 
^ ■ the framework of just the Standard Model, at which temperature did chemical 

h : 

equilibrium between quarks and gluons first become possible? 

In the present paper we address this question by considering carefully qq 
annihilation into gluons in the early universe. The lowest order process is qq — > gg, 
but one expects that when s ^> Aq CD , qq annihilation into many gluons becomes 
more and more important. Therefore in the early universe it might be essential 
to consider also qq — > gg ■ ■ ■ ^-annihilations. We show however that at very high 
temperatures these modify the annihilation rate by only about 3%. Comparing the 
total annihilation rate with the expansion rate of the universe, which we correct 
for the Standard Model interactions in the primeval plasma, we find that quarks 
and gluons were not in chemical contact above T > 3 x 10 GeV. We expect that 
at still higher temperatures the Standard Model particles did not have any thermal 



contact whatsoever, but precise study of kinetic equilibrium between e.g. gluons is 
made difficult by infrared problems and soft gluon physics, which cannot be dealt 
with succesfully in perturbation theory, and we do not address this issue here. 

Let us begin by computing the thermally averaged annihilation rate 
F&nniqq —> 99) 7 which in the limit when the final state blocking is neglected can 
be defined as 

T* m (qq^gg) = ±-f j0^j^^f(E 1 /T)f(E 2 /T)a( q q^gg)v rel p 1 .p 2 (1) 

where n q is the quark density (and equal to the anti-quark density), 
f(E/T) = (exp(E/T) + and v re \ = — (m 2 / (pi ■ P2)) 2 is the invariant 

velocity. The relevant t- and s-channel diagrams are depicted in Fig. 1. We find 
that 
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where x = s/4m 2 with m the quark mass, and A = CaTf(N 2 — 1) = N/2(N 2 — 1) 
and B = NC 2 F = 1/(4N)(N 2 - l) 2 are SU(N) color factors. Note that in the 
case of QED, A = and B = 1 so that Eq. (2) reduces to the well known Dirac 
formula. 

In the early universe it is natural to replace m by the quark plasma mass 
m q (T) in Eq. (1). At high temperature the left-handed quark and gluon plasma 
masses are given by [4] 

m%T) = 2 -g 2 T 2 , 

1 3 1 ^ 

ml{T) =(- 6 9>+ T2 9 2 W+ —gi Y )T 2 , 

where gw is the weak and gy the U(l)y coupling; for right-handed quarks the 
purely weak contribution is absent. Substituting these back to Eqs. (1) and (2) 
we find that T ann (qq — > gg) can be written as 

r„„( M -- S9 ) = ^|§^(iogr)T, (4) 



where F 2g (log T) is dimensionless quantity that can be solved numerically. Eq. (3) 
implies that F 2g can depend on T only through running of a s with temperature. 
Thus it is a smooth, slowly changing function of T and, neglecting SU(2) x U(l) 
corrections in the quark plasma masses, F 25 (3) = 431 and F 29 (15) = 680 (see 
Table 1). If all quarks were left-handed, we would obtain F^ 9 (3) = 397 and 
F 1 5 (15) = 597. In what follows, we do not compute the rate (4) separately for 
different chiralities but approximate F 2g by the average of F* 9 and Fl 9 . 

We should comment here that the QCD running coupling at finite temperature 
is, in general, a problematic concept, and there are several conflicting statements 
about the way it scales with temperature [5,6]. This has to do with the fact that 
at finite temperature, the system has two a priori independent mass scales, T and 
the renormalization point fi, so that the limit T — > oo is not unambigious. In our 
case thermal averaging, Eq. (1), tacitly assumes that the external legs are, apart 
from plasma mass corrections, free particles with momenta partitioned thermally. 
(This should be a good approximation for weakly interacting gas.) In that case 
the ensemble average G of the effective charge [1], defined as the thermal average 
of the usual effective charge, can be shown to scale as G(fx,T,g) — > G((j,,T,g) 
under T — > AT, where g = g(X/j,,T/ y) is the running coupling with T//jl fixed. 
Thus, in collisions the coupling is, on the average, the running finite T coupling [5] 
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with T//J, fixed. The origin of the function clq, which in the T//jl — > oo limit is 
a polynomial of T/fj,, are both the vacuum and the T-dependent finite parts 
appearing in the charge renormalization constant. That is, ao(T//x) gives the 
first order change to the running coupling g s when changing the renormalization 
scheme from the MS'-scheme: 

9 = 9ms (1 + a (T/iJ,)g 2 MS + 0(g 4 MS )) (6) 

Thus minimal sensitivity to the temperature corrections is obtained when ao 
vanishes. In the MOM-scheme for pure QCD this happens when \i ~ 2.6 T [1]. 
With this particular T//jl value guoM runs exactly like gMs{T = 0): 



In what follows we shall assume that one may also here adopt a prescription 
where a (T/n) = 0. Then we are able to convert the value of the strong coupling 
measured in LEP at ji ~ Mz and T = [7], a s = 0.12, to finite temperatures by 
using Eq. (5). In effect this means that g s (as well as gw and gy) runs exactly as 
in vacuum, but with fx replaced by kT with k a constant, and the LEP reference 
point is translated to the temperature T ~ Mz/k. 

The computation of qq — > ggg is more involved as compared with qq — > gg 
because of the more complicated kinematics of 2 <-> 3 scattering. The relevant 
diagrams and definitions of kinematical variables are shown in Fig. 2. We perform 
the calculation in large-N approximation where the leading part of the the matrix 
element for massless quarks and gluons, summed over colour and helicity degrees 
of freedom, can be written as [8] 

1 n 1 

\M(qq - ggg)\* = 2g 6 s N 2 (N 2 - 1)- £ (s%s qi + s qi s%) £ — - — - + 0(N~ 2 ) 

S ^ ^ S ql S 12 S 23 S 3q 

(8) 

where s^- = (pi + Pj) 2 , s q i = (p q — Pi) 2 and s q i = (p q — Pi] 2 ■ The second sum in 
Eq. (8) extends over permutations of the gluon indices. 

The 2^3 cross section can be written, in the notation of [9], as 

o^(jV 2 -l) f dt 1 dt 2 ds 1 ds 2 P^si, s 2 ,h,t 2 ) , ft * 
<r(qq -»• ggg) = ™ — 77 2?^ — 27^vT / ; : n tz > ( 9 ) 
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where A 4 is the Gram determinant defined in [9] and Sj and £j are invariants 
defined in Fig. 2. The term P 4 is the first sum in Eq. (8) with invariants 
written in terms of Sj and in the center-of-mass frame. Although the cross 
section Eq. (9) has infrared problems, in the early universe the plasma masses of 
quarks and gluons act as natural regulators. We account for them by replacing 
— m 2 in Eq. (9), where m is the relevant plasma mass. Note that 
plasmon decay into quarks is not kinematically possible so that 1 <-> 2 processes 
should not play a major role in equilibration. 



We have evaluated the cross section (9) numerically, but when s ^> m 2 , we find 
that it is well fitted by 
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The fit is valid when s/m > 10 . Comparing with a((/(/ — > gg) we find 
that annihilation into three gluons, taking m g = 2m g , is equally important 
when s/m 2 ~ 10 4 with a s ~ 0.12. For very large s/m 2 , 3-gluon final states 
begin slowly to dominate over 2-gluon final states. For example, when 
s/m 2 ~ 10 11 , a(qq -> ggg)/a(qq -> gg) ~ 5. 

In the early universe s in collisions was on the average about 18T 2 , not 
significantly larger than the squared thermal masses. Therefore one should not 
expect an enhancement of the average annihilation rate due to the large logarithms 
of the qq — > ggg cross section. We have evaluated the thermally averaged qq — > ggg 
rate as in Eq. (1), and we find numerically that 

r a „„(^ - ggg) = ^J—F s 9(\ogT)a 3 s (T)T (11) 

where F 3g is tabulated in Table 1, and we see that at very high temperatures 
2^3 rate is about 3% of the 2^2 rate. Thus, although quarks and antiquarks 
annihilate into any number of gluons in the early universe, these processes should 
not contribute significantly to the total annihilation rate T^ D . 

We now compare the total averaged annihilation rate T^ a with the expansion 
rate of the universe as given by the Hubble parameter H = ^/8tvGnp/3. In what 
follows we shall take into account the fact that the gas in the early universe 
is interacting, which modifies the energy density p from its ideal gas value. 
The thermodynamic potential O has been computed for SU(N) gauge theories in 
perturbation theory [10], and to lowest it reads 



7T 2 T 4 / 7 \ 

n = - — (N>-l + -NN f ), 

V J (12) 

1 144 I 4 ' 
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where fio is thermodynamic potential for ideal gas and fii is the lowest order 
exchange energy correction, and Nf is the number of flavours. For U(l) one 
should set N 2 - 1 -> 0, N -> 1 in fi and N 2 - 1 -> 1, A" = 0, A 7 / -> £ Qj in 
Energy density is given by p = fl + TS = 0- Tdfl/dT = -30 with O = fi + fii- 
In the Standard Model at T > 100 GeV we find that 



so that at high temperature psm differs from ideal gas value by a few percents. 
The effect of the exchange energy in the case of the Standard Model is illustrated 
in Fig. 3 for T ^> 100 GeV, assuming Higgs and top masses can be neglected. 

In Fig. 4 we have drawn H together with r^ n ~ T ann (qq -> gg)+T ann (qq -> ggg). 
We see that the total annihilation rate is too slow to maintain chemical equilibrium 
when 



The effects considered here work all in the same direction: including qq — > ggg 
in T)^ n increases the decoupling temperature Td, and so does also the exchange 
energy (12) by slowing down the Hubble rate. Annihilation into four or more 
gluons, which in principle could increase Td further, is unlikely to change the 
conclusion (14) because of the smallness of the logarithmic enhancement of the 
multiparticle production rate. Kinetic equilibrium, which is maintained by elastic 
collisions, may extend to temperatures somewhat higher than Td, but also kinetic 
contact is lost at sufficiently high temperature. Thus, unless the quark and 
gluon ensembles were, for some reason, created thermal at energy scales Mpi, 
they cannot have had thermal distributions at and above GUT scales. The 
same conclusion would naturally hold also for the rest of the Standard Model 
particles. This is directly relevant to the treatment of phase transitions at very 
high temperatures, which often is based on the assumption that the background 
Higgs field (or the scalar order parameter) evolves in a thermal background. In 
reality the familiar thermal corrections to the effective potential above Td may be 
completely absent, or present only in a modified form. 




(13) 



T = T d > 2.5 x 10 14 GeV. 



(14) 
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Table caption 



log T 


^0 


F 2 9 
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431 


397 
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6 


525 


474 
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590 


526 


10 


12 


640 


565 


15 


15 


680 


597 


19 



Table 1. The behaviour of the functions F 9 , F x 9 and F 39 with temperature. 
Fq 9 contains only the SU(3) C contribution to the thermal particle masses, while 
F^ 9 accounts also for the electroweak corrections but assumes that all quarks are 
left-handed. F^ 9 has only SU(3) C corrections. 
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Figure captions 



Figure 1. The annihilation diagrams for qq — ► gg. 

Figure 2. The annihilation diagrams for qq — > ggg, and the definitions of the 
relevant kinematical variables. 

Figure 3. The Standard Model degrees of freedom at high temperatures, corrected 
for the lowest order exchange energy (solid curve). For comparision, we show also 
the ideal gas value g* = 106.75 (dashed line). 

Figure 4. Comparision of the Hubble rate (solid curve), as computed from (13), 
with the qq annhilation rate T^ n (dashed line). 
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